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Abstract 

The National Aeronautics and Space Administration (NASA) is currently investigating a 
conceptual robotic mission to collect a small boulder up to 4 m in diameter resting on the 
surface of a large Near Earth Asteroid (NEA). Because most NEAs are not well characterized, 
a great range of uncertainties in boulder mass properties and NEA surface characteristics must 
be considered in the design of this mission. These uncertainties are especially significant when 
the spacecraft ascends with the boulder in tow. The most important requirement during ascent 
is to keep the spacecraft in an upright posture to maintain healthy ground clearances for the 
two large solar arrays. This paper focuses on the initial stage (the first 50 m) of ascent from 
the surface. Specifically, it presents a sensitivity study of the solar array ground clearance, 
control authority, and accelerations at the array tips in the presence of a variety of uncertainties 
including various boulder sizes, densities, shapes and orientations, locations of the true center 
of mass, and push-off force distributions. Results are presented, and appropriate operations are 
recommended in the event some of the off-nominal cases occur. 


1 Introduction 

In recent years the National Aeronautics and Space Administration (NASA) has been considering 
a conceptual robotic mission to collect a small Near-Earth Asteroid (NEA) or a small boulder 
resting on the surface of a large NEA, and transport it to an orbit in the Earth-Moon system. The 
latter option is referred to as the Asteroid Robotic Redirect Mission (ARRM) Alternate Concept 
[1]. Such a mission would lead to greater scientific understanding of NEAs and could demonstrate 
concepts for planetary defense. Dynamics and control during various phases of the mission when 
the spacecraft is in the vicinity of the asteroid are discussed in Refs. [2] and [3]. Specifically, 
initial approach, hover, spacecraft spin-up, descent, ascent (in the case of boulder retrieval), and 
de-spin are examined in order to determine the velocity increments and control force and torque 
that must be provided by a Reaction Control System (RCS), and the mass of the propellant that 
will be consumed. The present paper is concerned with retrieving a boulder from the surface of 
a large NEA; in particular, the focus is on ascent from the NEA and the dynamics and control of 
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the trajectory and attitude of the spacecraft carrying a boulder that has been secured by a capture 
mechanism. 

As discussed in Ref. [4], the current design of the ARRM Alternate Concept includes three 
mechanical legs to contact the surface and stabilize the spacecraft, and provide the liftoff forces 
necessary to escape the NEA gravity field. In addition, three manipulator arms are used to secure 
the boulder; in this study we assume that the spacecraft, manipulator arms, and boulder form one 
composite rigid body during ascent. After the push-off legs leave the surface, RCS jets are used 
to control the attitude. The capture mechanism is designed to handle a baseline boulder diameter 
of 4 m, and can accommodate a diameter as small as 0.5 m. The primary goal of the controlled 
ascent is to depart from the NEA and maintain attitude control throughout the first 10 to 50 m such 
that the solar array tips maintain a healthy ground clearance. The wingspan of the spacecraft is 
approximately 35 m with the MegaFlex arrays and 51m with the Rosa Arrays. In this work, the 
MegaFlex arrays are considered. 

For the purpose of this study, the NEA target is selected to be asteroid 25143 Itokawa, which 
spins at 12.1324 hrs per revolution, and the landing site is chosen to be in the area of the Muses 
Sea (Ref. [5]). A three-axis attitude controller is designed to track the local vertical in real time, 
with the goal of keeping the solar array axes parallel to the ground during the early stage of ascent. 
This controller is based on the inertia- free linear feedback controller presented in Ref. [6], which 
is robust to the unknown moment of inertia of the rigid body if the control torque is provided by an 
actuator that exerts a pure couple, such as a reaction wheel or other momentum exchange device. 

The absence of good physical characterizations of the majority of NEAs and their boulders 
makes controlled ascent a challenging prospect. Little time is available for boulder mass property 
estimation during the early portion of ascent because it is important to leave the vicinity of the NEA 
quickly in order to avoid a collision. This paper presents sensitivity analysis of attitude control 
performance in the presence of uncertainties in the following parameters: boulder size, estimated 
density, density uncertainty, shapes and orientations of the boulder relative to the spacecraft when 
captured, the unknown offset of the true center of mass (CM) from the estimated location, and 
uneven distribution of force among the three push-off legs. The metrics of primary interest are the 
deviation of the array axis from the local horizontal plane, the ground clearance of the solar array 
tips, the acceleration experienced at the solar arrays, and the control authority. The results of a 
large number of numerical simulations show that attitude control performance is most affected by 
CM offset and an uneven distribution of push-off forces. 

To the best of the authors’ knowledge, controlled ascent has not been performed from a celes- 
tial body other than the Moon except for an emergency ascent of the Hayabusa spacecraft from 
Itokawa, which became known only later when data was analyzed (Ref. [5]). No precedent exists 
for controlled ascent with a payload whose size and mass are comparable to those of the spacecraft, 
and whose physical characteristics are not well known. 

The remainder of this paper is organized as follows. A description of the spacecraft, landing 
legs, and capture arms is provided in Sec. 2, and an ascent Concept of Operations that is particular 
to the study is presented in Sec. 3. The sources of uncertainties and the trade space are discussed 
in Sec. 4, the main results are presented in Sec. 5, and conclusions are given in Sec. 6. 
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2 The Spacecraft 


The spacecraft is shown with a boulder in Figs. 1 and 2. The spacecraft bus is roughly a right 
hexagonal prism with a mass of 7,600 kg, a height of 4.96 m, and a diameter of 3 m. Two solar 
arrays are attached to the bus; each has a mass of 200 kg, and a center located 11.805 m from 
the centerline of the bus. The spacecraft lands on the asteroid with three legs, which support the 
spacecraft during its stay on the asteroid’s surface. When the spacecraft lands on the asteroid’s 
surface, the legs pads are 3 m away from the centerline. Each leg consists of multiple articulated 
truss-structure links. Three similar mechanisms, smaller than the legs, are used to collect the 
boulder and secure it to the spacecraft prior to ascent and departure from the asteroid’s surface. 
When the boulder is captured, it is held such that the boulder CM is 3.5 m away from the base of 
the capture arms. A comprehensive description of the landing legs and capture arms can be found 
in Ref. [4] . 

Four pods, each containing four Reaction Control System (RCS) thrusters, are mounted to the 
bus on struts. The position vector r B p ‘ from the mass center of the bus, B*, to each of the four pods 
Pi (i = 1,2, 3, 4), and direction of thrust applied by each thruster, are reported in Table 1 in terms 
of bi, 62, and 63, a set of right-handed mutually orthogonal unit vectors fixed in the spacecraft 
bus, directed as shown in Figs. 1 and 2. Each thruster has a specific impulse 7 sp = 230 s. The force 
applied by thrusters 1, ... ,4 has a maximum magnitude of 15.6 N, whereas the force applied by 
thrusters 5, . . . , 16 has a maximum magnitude of 22.2 N. 
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Figure 1 : Side view of spacecraft with boulder 



Figure 2: Bottom view of spacecraft with boulder 
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Table 1 : RCS Thruster Configuration 


Pod 


Position (m) 

No. 

Thrust Directions 

Pi 

r b*P\ _ 

1.5bi + 1.5b 2 -2.48b 3 

1 

b 3 




5 

-b 2 




9 

1 /v 1 /V 

->-72 b * 




13 

1 /v 1 /V 

- 72 ^- 72 ^ 

Pi 

r B *Pi = 

— 1.5bi + 1.5b 2 — 2.48b 3 

2 

b 3 




6 

— b 2 




10 





14 

-72^-72*1 

Pi 

r FPi = 

— 1.5bi — 1.5b 2 — 2.48b 3 

3 

b 3 




7 

b 2 




11 

72^ + 72^ 




15 

^b 2 -^b 3 

Pa 

r B * P 4 - 

1.5bi — 1.5b 2 — 2.48b 3 

4 

b 3 




8 

b 2 




12 

1 /v 1 /V 

-> + /5 b 2 




16 

1 /v 1 /V 

72^-72^ 


3 Ascent Concept of Operations 

The proximity operations of the ARRM Alternate Concept consists of asteroid characterization, 
dry run attempts, approach, descent, landing, capture, ascent, etc., which are well described in Ref. 
[7]. For this study of the ascent phase, it is assumed that the landing site is located in the Muses Sea 
at coordinates (126.3, 76.8, —12.2) m in the body frame of Itokawa [8]. The axis of symmetry of 
the spacecraft is assumed to be initially aligned with the local vertical. The push-off legs provide a 
constant force to achieve a AF = 0.2 m/s in 10 seconds, which is enough to escape the gravity field 
of Itokawa. Due to possible uneven surface firmness beneath the three leg pads, or uncertainties 
in the actuators, the legs may provide uneven push-off forces. The force from each leg is assumed 
to remain constant throughout the 10 seconds. The uncertainties in boulder mass properties and 
CM position, together with uneven push-off forces, can lead to a disturbance torque that causes 
the solar arrays tip away from local horizontal and dip towards the surface. This can happen at 
the very beginning of the ascent; therefore, the attitude controller is engaged at all times, with the 
commanded torque supplied by the RCS thrusters. In addition, it is assumed that the cohesive force 
that keeps the boulder attached to the asteroid is already broken prior to this part of the analysis, 
and the dynamics of liberating the boulder is not considered. 

In order to avoid pluming the NEA surface and stirring up dust that may contaminate the 
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spacecraft, thrusters 13, 14, 15 and 16 are not used during the first 20 m of ascent. In addition, 
mission requirements state that the ground clearance of the solar arrays must be above 5 m, and 
the acceleration at the tips of the arrays must be less than 0.1 g. 


4 The Sources of Uncertainties and Trade Space 

Major challenges to the success of the mission are presented by a number of uncertainties during 
the initial stage of ascent, including the mass properties of the boulder, the unknown characteristics 
of the landing surface, navigation errors, unevenly distributed push-off forces, and CM offset that 
can be caused by a number of factors. All of these parameters greatly affect the performance 
of the controller designed to maintain an upright posture of the composite body, and effects of 
parameter variations are studied in this paper. The parameters of interest are boulder size, boulder 
density, boulder density uncertainty, boulder shape, boulder orientation relative to the array axis, 
CM offsets, and push-off force distribution. The total number of cases in the trade space is 8100. 
The uncertainties assumed in this study are now described and quantified. 

Several combinations of boulder size and shape are studied. Three shapes are considered, 
namely a sphere, a 21 1 ellipsoid, and a 221 ellipsoid. Ratios of diameter lengths are indicated in the 
numerical designations of the ellipsoids. For each shape, three sizes are considered, corresponding 
to a spherical boulder diameter d of 4 m, 3 m, and 2 m. The diameters of the two ellipsoids are 
determined by assuming that the ellipsoid has the same mass as the sphere with a given diameter. 
The semi-diameters of the 211 and 221 ellipsoids are listed in Table 2. 


Table 2: Semi-diameters a, b, c, (m) of the considered boulders 


Size 

Sphere 

211 Ellipsoid 

221 Ellipsoid 

a 

b 

c 

a 

b 

c 

a 

b 

c 

d — 4 m 

2 

2 

2 

3.17 

1.59 

1.59 

2.52 

2.52 

1.26 

d — 3 m 

1.5 

1.5 

1.5 

2.38 

1.19 

1.19 

1.89 

1.89 

0.94 

d — 2 m 

1 

1 

1 

1.59 

0.79 

0.79 

1.26 

1.26 

0.63 


In the case of the 221 ellipsoid, it is assumed that the boulder sits on the surface with its short 
axis perpendicular to the surface. When the boulder is a 21 1 ellipsoid, it is assumed that the boulder 
lies on the surface with its long axis parallel to the surface. For the 211 ellipsoidal boulder, two 
orientations of the long axis are considered: parallel and perpendicular to the solar array axis. 

It is assumed that a good shape model of the boulder can be established prior to the capture 
during the characterization phase and dry runs to the candidate target sites. It is also assumed 
that the boulder is coherent and that there is no significant variation in the density throughout 
the boulder. Based on these assumptions, the center of mass and the moment of inertia of the 
boulder can be calculated, provided the density of the boulder can be estimated. In this study, it 
is assumed that the boulder density can be estimated to within ±10% of the true density. Three 
values of density are considered, namely, 3.0, 2.5, and 2.0 g/cc. For each value, the true density 
is considered to be 10% less than, the same as, or 10% more than the estimated density. The 
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Table 3: Composite body vertical CM offset (mm) along $3 axis (relative to the estimated position) 
due to density variation 


Est. density 

2.0 g/cc 

2.5 g/cc 

3.0 g/cc 

True density 

-10% 

+ 10% 

-10% 

+ 10% 

-10% 

+ 10% 

d = 4 m 

-6.3 

5.2 

-5.2 

4.4 

-4.5 

3.7 

d = 3 m 

-11.1 

9.5 

-9.8 

8.3 

-8.7 

7.4 

d = 2 m 

-15.7 

14.2 

-15.6 

13.9 

-15.1 

13.4 


Table 4: Boulder masses (ton) considered in the trade space 


Diameter 

2.0 g/cc 

2.5 g/cc 

3.0 g/cc 

s/c mass 

d = 4 m 

67.0 

83.8 

100.5 


d = 3 m 

28.3 

35.3 

42.4 

8 

d = 2 m 

8.4 

10.5 

12.6 



uncertainties in the density estimation affect not only the moment of inertia of the boulder, but also 
the mass of the boulder, which consequently causes uncertainties in the location of the CM of the 
composite body. 

For spherical boulders, the effects of variations in diameter and density on the offset of the 
combined CM (true CM relative to the estimated CM) in the direction of 63 are recorded in Table 
3. It can be seen that the vertical CM offset increases when the boulder gets smaller and lighter, 
and up to 16 mm of vertical CM offset can be expected for the smallest boulder considered. The 
CM offset affects the knowledge of the combined moment of inertia and the accuracy of the torque 
delivered by the RCS thrusters. 

With the boulder sizes and densities, the boulder mass ranges from being comparable to the 
spacecraft mass to up to two orders of magnitude heavier. The boulder mass for the cases consid- 
ered in this study as well as the spacecraft mass are shown in Table 4. 

The spacecraft uses three supporting legs to provide the push-off force needed to achieve an 
escape velocity of 0.2 m/s off the surface of Itokawa. However, unknown surface characteristics as 
well as uncertainties in the push-off actuators themselves can contribute to uneven distribution of 
push-off forces among the three legs. Five force distributions considered in this study are listed in 
Table 5. 

When the spacecraft touches down on the surface, its axis of symmetry may not be accurately 
aligned with the CM of the boulder due to navigation errors during the descent. The position error 
of the onboard sensors is taken to be 0.5 m, 3a. This contributes to the knowledge errors of the 
true CM location of the composite body. In this study, the sensitivity of the ascent performance 
with respect to the CM offset is investigated, and the trade space includes boulder CM offsets of 
0.5 m in the 15 directions listed in Table 6. 
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Table 5: Push-off force distributions among three legs 


Index 

leg 1 

leg 2 

leg 3 

1 

Even distribution 

2 

30% 

40% 

30% 

3 

40% 

30% 

30% 

4 

25% 

50% 

25% 

5 

50% 

25% 

25% 


Table 6: Location of true boulder CM relative to estimated CM in the trade space (unit = m) 


Index 

CM offsets 

1 

no offset 

2 


3 

4 b i 

4 

\b 2 

5 

~\b 2 

6 

2 b 3 

7 

4 b i 

8 

#bi + 4b 2 + 4 b 3 

9 

#b,-f6 2 + fb 3 

10 

-^b,-fb, + fb 3 

11 

-#bi + 4b 2 + 4 b 3 

12 

#b, + fb 2 -#b 3 

13 

^b,-fb 2 -fb 3 

14 

-#b,~fb 2 -fb 3 

15 

-^bi + 4 b 2-4 b 3 
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5 Results 


This section summarizes the results of the sensitivity trade study. A polyhedral gravity model [9] 
based on a shape model consisting of of 25,350 vertices and 49,152 facets is used for Itokawa. 
Since the force supplied by the push-off legs is much greater than that due to the solar radiation 
pressure, the effect of the solar radiation is neglected. 

A three-axis attitude controller is used to align the spacecraft axis of symmetry with the local 
vertical. It is based on the linear quaternion feedback controller presented in Ref. [6]. This ap- 
proach is chosen because the controller does not rely on the knowledge of the moments of inertia 
of the rigid body, and it is thus an inertia-free feedback control law. If the spacecraft is equipped 
with a torque actuator (an actuator that applies a pure couple), then the controller is robust to 
the uncertainties in the mass properties of the rigid body, provided enough control authority is 
available. In the present case the control torque is provided by the RCS thrusters and appropriate 
thrusters are selected to deliver the commanded torque. However, the thruster selection is based on 
the estimated location of the CM of the combined body. Thus, the actual torque delivered by the 
selected thrusters will be different from the commanded torque because of the uncertainties in the 
knowledge of the system CM. 

The knowledge of the spacecraft attitude and angular rate is provided by the onboard navigation 
system, which is not perfect. In this study, the standard deviation for the attitude error is assumed 
to be 0.01 deg (36 arcsec), and the standard deviation for the angular rate error is assumed to be 
1/100 of the Itokawa spin rate (1.4 x 10 6 rad/s). 

5.1 The boulder size and density 

The effect of the uncertainties in the trade space greatly varies with the size and density of the 
boulder. The larger and heavier the boulder is, the harder it is in general to control the composite 
body in the presence of perturbations. For the boulder diameters and densities considered, the 
maximum angles between the solar array axis and and local horizontal plane are listed in Table 
7, and the minimum ground clearances of the solar array tips are shown in Table 8. As can be 
seen in the tables, as the size or density of the boulder increases, the solar array experiences larger 
deviations from the local horizontal plane, resulting in smaller ground clearances. 


Table 7: Maximum deviation (deg) of the array axis from the local horizonal plane 


Diameter 

density 

2 g/cc 

2.5 g/cc 

3 g/cc 

2 m 

3.05 

3.82 

4.73 

3 m 

15.74 

23.07 

30.49 

4 m 

66.56 

89.94 

89.72 


For a boulder of 4 m in diameter and all three values of density, in the worst case of trade space 
the solar arrays hit the ground while ascending as indicated by the negative ground clearances in 
Table 8. Table 7 shows that the spacecraft can tip over so much that the solar array axis is almost 
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Table 8: Minimum ground clearance (m) of the array tips 


Diameter 

density 

2 g/cc 

2.5 g/cc 

3 g/cc 

2 m 

9.52 

9.51 

9.50 

3 m 

8.29 

6.94 

5.13 

4 m 

-1.91 

-6.41 

-9.05 


perpendicular to the surface. The situation improves greatly when the boulder diameter reduces 
to 3 m. In this case, even though the solar arrays deviates from the local horizontal plane up 
to 30.49 degrees in the worst case, the ground clearance of the array tip is above 5 m, which is 
acceptable in terms of performance. It should be noted here that as the diameter, and hence the 
mass, of the boulder increases, the force per leg increases to achieve the 0.2-m/s escape velocity in 
10 seconds. Thus, the total disturbance torque applied to the system as a result of uneven push-off 
force distribution grows as the boulder becomes bigger and heavier, while the RCS thrusters are 
kept at the same level. 



Figure 3: The deviation of the solar array 
axis from the local horizontal for the case 
with the worst ground clearance with a 4-m 
diameter boulder. 



Figure 4: The worst-case ground clearance 
of the solar arrays with a 4-m diameter boul- 
der. 


For the case of a 4-m diameter boulder, the deviation of the array axis from the local horizontal, 
ground clearance of the arrays, and the control torque for the case with the worst ground clearance 
are shown in Figs. 3-5. With a AV — 0.2 m/s achieved within 10 seconds, it takes about 227.5 
seconds for the combined body to ascend to a 50-m altitude. Figure 3 indicates that the solar 
arrays tip past the local vertical, and does not level out throughout the 50-m ascent. Figure 4 shows 
that one of the arrays hits the surface when the combined body is at an altitude of 7.5 m, about 
37.5 s after the start of the ascent. Figure 5 shows that in this case, the commanded torque is 
significantly greater than the torque that can be delivered by the RCS thrusters, causing the system 
to be severely underactuated. 
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Figure 5: The control torque profile for the case with the worst ground clearance with a 4-m 
diameter boulder. 

For the case of a 3-m diameter boulder, the deviation of the array axis from the local horizontal, 
ground clearance of the arrays, and the control torque for the case with the worst ground clearance 
are shown in Figs. 6-8. Figure 6 indicates that the solar arrays are leveled out about 150 s into 
the ascent, at which time an altitude of approximately 30 m is reached. During the ascent, the 
spacecraft tips over while leaving the surface, so the deviation from the local horizontal plane is 
not the only indicator of the array ground clearance. In fact, Fig. 6 indicates that the maximum 
deviation angle occurs around 40 s, or 8 m altitude, but Fig. 7 shows that the arrays are closest to 
the ground at about 5 m in altitude, or about 25 s after the start of the ascent. Figure 8 shows that 
in this case, the commanded torque is significantly greater than the torque that can be delivered 
by the RCS thrusters. Nonetheless, the RCS thrusters have enough control authority to keep the 
arrays off the ground by at least 5 m in this case. 

If the boulder is 2 m in diameter, then the worst case deviation of the solar array axis from the 
local horizontal plane is less than 5 deg, with a minimum ground clearance of 9.5 m, as shown in 
Tables 7 and 8. Notice that the arrays have a ground clearance of 9.52 m when the spacecraft is on 
the ground. Thus, in the worst case, the arrays only dip 20 mm below the initial level, which leaves 
a safety margin of 4.5 m. For the case of a 2-m diameter boulder, the deviation of the array axis 
from the local horizontal, ground clearance of the arrays, and the control torque corresponding to 
the case with the worst array ground clearance are shown in Figs. 9-11. Figure 9 indicates that the 
solar arrays are essentially leveled out about 60 s into the ascent, at which time an altitude of about 
12 m is reached. Figure 11 shows that in this case, the commanded torque can be 2.5 times greater 
than the torque that can be delivered by the RCS thrusters. Nonetheless, the RCS thrusters have 
enough control authority to keep the arrays off the ground by 5 m with a healthy margin. 
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Figure 6: The deviation of the solar array Figure 7: The worst-case ground clearance 
axis from the local horizontal for the case of the solar arrays with a 3-m diameter boul- 
with the worst ground clearance with a 3-m der. 
diameter boulder. 



Figure 8: The control torque profile for the case with the worst ground clearance with a 3-m 
diameter boulder. 

5.2 Boulder shape and orientation 

Boulder shape affects the moments of inertia of the combined body, which significantly influences 
the attitude control performance. Table 9 lists the principal moments of inertia of 4-m diameter 
boulders with a density of 3 g/cc in the shapes considered. It can be seen that In and I 22 for the 
spherical and 221 ellipsoidal shapes are virtually the same (less than 0.8% different), and only In 
is significantly different. Since it is predominantly 1 1 1 and I 22 that affect the control of the array 
axis orientation, it can be expected that the array ground clearances do not vary much when a 
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Figure 9: The deviation of the solar array 
axis from the local horizontal for the case 
with the worst ground clearance with a 2-m 
diameter boulder. 



Figure 10: The worst-case ground clearance 
of the solar arrays with a 2-m diameter boul- 
der. 



Figure 1 1 : The control torque profile for the case with the worst ground clearance with a 2-m 
diameter boulder. 

boulder assumes a spherical shape or a 221 ellipsoidal shape. In fact, the results show that array 
ground clearances for a spherical boulder and a 221 ellipsoidal boulder differ by no more than 20 
mm. In addition, in the case of a spherical boulder and a 221 ellipsoidal boulder with the short 
axis perpendicular to the surface, the relative orientation about the axis of symmetry between the 
spacecraft and boulder is irrelevant. 

The principal moments of inertia are significantly different between a sphere and an equivalent 
211 ellipsoid with the same mass, resulting in vastly different performances. Figure 12 shows the 
comparison of the array ground clearances for a spherical boulder, a 21 1 ellipsoidal boulder with 
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its long axis parallel to the array axis, and a 21 1 ellipsoidal boulder with its long axis perpendicular 
to the array axis. For the cases shown, the boulder diameter is 4 m, the estimated density is 3 g/cc, 
and the true density is 3.3 g/cc. It can be seen from the figure that for all combinations of CM offset 
and push-off force distributions, the case of the 211 ellipsoidal boulder with its long axis parallel 
to the array axis has the greatest array ground clearances, followed by the case of the spherical 
boulder, and the case of the 211 ellipsoidal boulder with its long axis perpendicular to the array 
axis has the smallest array ground clearances. 

This result can be explained by comparing the moments of inertia of the combined body about 
the 62 axis, which is perpendicular to the array axis. It is clear that the combined moment of inertia 
about the 62 axis is the smallest for a 211 ellipsoidal boulder with its long axis perpendicular to 
the array axis, and the greatest for the same boulder when its long axis is parallel to the array axis. 
A greater moment of inertia about 62 axis results in a slower response given the same disturbance 
torque, and consequently greater array ground clearance compared to a smaller moment of inertia 
about the t >2 axis. 


Table 9: Principal moments of inertia (kg-m 2 ) of 4-m diameter boulders with density of 3 g/cc of 
various shapes 


Shapes 

In 

I22 

I33 

Sphere 

1.61 x 10 5 

1.61 x 10 5 

1.61 x 10 5 

221 ellipsoid 

1.60 x 10 5 

1.60 x 10 5 

2.55 x 10 5 

211 ellipsoid 

1.01 x 10 5 

2.53 x 10 5 

2.53 x 10 5 


5.3 Density uncertainty 

Uncertainty in density estimation determines whether the real boulder is heavier or lighter than the 
estimated value. It affects not only the accuracy of the estimate of the moment of inertia of the 
composite body, but also the location of the combined CM, which in turn affects the accuracy of the 
torque supplied by the RCS thrusters. Figure 13 shows the array ground clearances for different 
density estimate uncertainties and combinations of CM offsets and push-off force distributions. 
This is the case corresponding to a 4-m 211 ellipsoidal boulder with its long axis perpendicular 
to the array axis and an estimated density of 3 g/cc. It can be seen that for all CM offsets and 
push-off force distributions, a true density greater than the estimated density results in a decrease 
of the ground clearance of the arrays, and the opposite is true if the true density is smaller than 
the estimated density. This trend holds true for all other combinations of boulder size, shape, 
orientation, and density. 

The maximum differences in array ground clearances between the cases with density uncertain- 
ties and the cases with no density uncertainties for the five assumed push-off force distributions 
are shown in Table 10. It can be seen that the effect of density uncertainty is most prominent in 
the cases of more extreme uneven push-off distribution (50%-25%-25%), and less so in the even 
push-off cases. A density uncertainty of 10% can result in up to 2 m of ground clearance difference 
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Figure 12: The array ground clearances for a 4-m diameter boulder with density 3 g/cc in different 
shapes and orientations, (a) even push-off, (b) 30%-40%-30%, (c) 25%-50%-25%, (d) 50%-25%- 
25%, (e) 40%-30%-30%. 

in the presence of uneven push-off force distribution. This data clearly shows that overestimating 
the density results in greater array ground clearance than underestimating the density. 

5.4 CM offset and uneven push-off force distribution 

The greatest contribution to the off-nominal performance of the spacecraft ascent is undoubtedly 
from the combination of the CM offset and the uneven push-off force distribution. CM offset con- 
tributes to the uncertainty in the estimation of the true moment of inertia and the torques delivered 
by the RCS thrusters, while the uneven push-off force distribution contributes to the attitude dis- 
turbance during the push-off. The worst-case array ground clearances for the combinations of CM 
offsets and push-off force distributions are shown in Fig. 14. This corresponds to the case involv- 
ing a 4-m 211 ellipsoidal boulder with its long axis perpendicular to the array axis, an estimated 
density of 3 g/cc, and a true density of 3.3 g/cc. 

It can be seen that in the case of an even push-off, the minimum array ground clearance is 

8.4 m among all CM offsets, with a healthy 3.4 m margin above the 5-m safety array clearance 
requirement. On the other hand, in the event of no CM offset (CM offset index = 1), all uneven 
push-off force distributions except 50%-25%-25% provide at least a 3.7 m margin above the 5-m 
clearance requirement. In the case that is the exception, the array clearance is 1.9 m and below the 
requirement of 5 m. 
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(d) 


(e) 




Figure 13: The array ground clearances for a 4-m diameter boulder with density 3 g/cc with dif- 
ferent density uncertainties, (a) even push-off, (b) 30%-40%-30%, (c) 25 %- 50 %- 25 %, (d) 50%- 
25%-25%, (e) 40%-30%-30%. 

It is observed from Fig. 14 that for push-off force distribution of 30%-40%-30%, ground clear- 
ances for all cases of CM offsets are above 5 m, with the worst case array ground clearance at 5.82 
m, a healthy margin of 0.82 m above the requirement. When leg 1 supplies 40% of the push-off 
forces, as in the case of 40%-30%-30%, the arrays on all but one occasion have ground clearance 
above 5 m with at least a 1-m margin. The exceptional case has a ground clearance of 1.7 m. 

Unsurprisingly, worse ground clearances are observed when the push-off forces are distributed 
more unevenly. If the push-off force is distributed 25%-50%-25% between the three legs, the 
array tips dip below 5 m to the surface in four occasions of CM offsets. When the push-off force 
distribution becomes 50%-25%-25%, all but five CM offsets result in array ground clearances less 
than 5 m. In five occasions, the arrays hit the ground. It should be noted that the worst ground 
clearance is observed when the true CM is located at — 0.5b i, and the leg 1 supplies 50% of 
the push-off force. In this case the true CM and the leg that supplies the largest push-off force are 
located on opposite sides of the spacecraft axis of symmetry, resulting in the maximum disturbance 
torque during push-off. 

The worst case in the trade space in terms of array ground clearance consists of trade parameters 
listed in Table 11. For this case, the accelerations experienced at the solar array tips are shown in 
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Table 10: Maximum difference in ground clearance (m) between the cases with uncertain densities 
and the cases with nominal density 


True 

Push-off force distribution (%) 

density 

even push-off 

30-40-30 

25-50-25 

50-25-25 

40-30-30 

-10% 

0.20 

0.65 

1.75 

1.98 

1.22 

+ 10% 

-0.25 

-0.75 

-1.76 

-1.86 

-1.32 



Figure 14: Worst case array ground clearance for combinations of CM offsets and push-off force 
distributions. 

Fig. 15. It can be seen that the largest solar array tip acceleration is less than 0.007 g, well within 
the 0.1 -g threshold. It should be noted that this analysis assumes a rigid array structure, and thus, 
these accelerations are likely minimum values that will be experienced. The large flexible structure 
of the arrays will encounter some whipping effect that will increase the peak values. 

5.5 Thruster sizing 

Analysis is also performed investigating how the current RCS thrust level has to be scaled in order 
to achieve array ground clearance of 5 m for the entire trade space. The required scale factors for 
the five push-off force distributions are listed in Table 12. It can be seen that in the worst case 
when the push-off force distribution is 50%-25%-25%, 4.1 times the current RCS thrust level is 
required to achieve a 5-m array ground clearance. If the push-off force distribution is 25%-50%- 
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Table 1 1 : Trade parameters in the worst case 


Parameter 

Value 

Boulder diameter 

4 m 

Estimated density 

3 g/cc 

True density 

3.3 g/cc 

Boulder shape 

211 ellipsoid 

Boulder orientation 

Long axis perpendicular to array axis 

CM Offset 

— 0.5bi 

Push-off force distribution 

50%-25%-25% 



100 150 

Time (s) 


Figure 15: Worst case acceleration experienced by the tips of the solar arrays. 

25%, a scale factor of 2.0 is required, and if the push-off force distribution is 40%-30%-30%, a 
scale factor of 1.5 is required. 


6 Conclusions 

The ascent of a spacecraft retrieving a boulder up to 4 m in diameter from the surface of an asteroid 
is examined with regards to sensitivity of solar array ground clearance, attitude control authority, 
and accelerations at the array tips, to uncertainties in several parameters. Various boulder sizes, 
densities, shapes and orientations, CM offsets, and push-off force distributions are considered. It 
is shown through numerical simulation that with a 2-m diameter boulder, the spacecraft attitude 
can be controlled such that the solar arrays have a ground clearance greater than 9.5 m, with a 
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Table 12: Thruster scaling to keep array ground clearance above 5 m 



Push-off force distribution (%) 


even push-off 

30-40-30 

25-50-25 

50-25-25 

40-30-30 

Scale factor 

lx 

lx 

2. Ox 

4.1 x 

1.5 x 


4.5-m margin over the 5-m safety requirement. For a 3-m diameter boulder, the solar array ground 
clearances are also above the 5-m requirement, but the margin is reduced to 0.13 m. For a 4-m 
diameter boulder, there are cases in which the solar array ground clearances are less than 5 m, 
and cases in which the arrays hit the ground. It is found that overestimating the density results in 
greater array ground clearance than underestimating the density. The cases involving a spherical 
boulder and a 221 ellipsoidal boulder have similar array ground clearances with differences less 
than 20 mm. For a 211 ellipsoidal boulder, orienting the boulder’s long axis parallel to the array 
axis results in a larger array ground clearance than a perpendicular orientation. In fact, the largest 
array ground clearance is obtained with a 211 ellipsoidal boulder in the parallel orientation. If 
the push-off force is distributed evenly or as 30%-40%-30%, then all cases result in array ground 
clearances greater than 5 m. It is found that the array tip accelerations are less than 0.007 g, 
well within the 0.1 -g requirement. In order to keep the array ground clearance above 5 m, the 
RCS thrust level needs to be scaled by 1.5, 2.0, and 4.1 times for push-off force distributions of 
40%-30%-30%, 25%-50%-25%, and 50%-25%-25%, respectively. 
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